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On the Solutions of Linear Non-homogeneous Partial 

Differential Equations, 

By L. L. Steimxey. 



The integrals of a partial differential equation of the first order were first 
classified by Lagrange, who separated them into three groups; namely, the 
general, the complete, and the singular integrals. For a long time this classi- 
fication was thought to be complete. In fact, Forsyth, in his "Differential 
Equations/' published first in 1885, gives a supposed proof of a theorem stating 
that every solution of such a differential equation is included in one or other 
of the three classes named. This error is also carried through the second and 
third English editions and the two German editions, the last one being pub- 
lished in 1912. 

In 1891 Goursat pointed out, in his " Equations aux derivees partielles du 
premier ordre," that solutions exist which do not belong to any of these three 
classes, and showed, indeed, that the existing theory is not complete even for 
the simplest forms. 

In November, 1906, Forsyth, in his presidential address to the London 
Mathematical Society, emphasized the fact that the theory is incomplete, and 
in his closing remark said : " It appears to me that there is a very definite 
need for a re-examination and a revision of the accepted classification of inte- 
grals of equations even of the first order: in the usual establishment of the 
familiar results, too much attention is paid to unspecified form, and too little 
attention is paid to organic character, alike of the equations and of the inte- 
grals. Also, it appears to me possible that, at least for some classes of equa- 
tions, these special integrals may emerge as degenerate forms of some semi- 
general kinds of integrals ; but it is even more important that methods should 
be devised for the discovery of these elusive special integrals." Forsyth also, in 
an address delivered at the Fourth International Congress of Mathematicians, 
takes advantage of the opportunity offered, to again emphasize the incomplete- 
ness of the existing theory of partial differential equations of the first order. 

Some indication of the nature of the incompleteness can be seen from the 
equation 
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A general solution by Lagrange's method is afforded by 

/(*!,**) =0, 

z z 

where / is an arbitrary function and where <p!= — and <p 2 — • We can easily 

x y 

verify that 

y 

affords a solution of the original equation; but this solution cannot be put in 
the form of the Lagrange general integral, since the Jacobian 

d(z, x, y) 

is not identically zero. Consequently the so-called general solution is not a 
general solution in the true sense of the term, since not every solution of the 
differential equation can be put in that form. 

In the present paper we shall deal with the linear non-homogeneous equation 

ix iP ,= Z, p ( =*., (1) 

where Xi and Z are functions of the variables z, x x , x 2 , . . . ., x n . .We shall 
assume, as we may without loss of generality, that all common factors of Z, 
X, , X 2 ,...., X n have been removed ; and consequently we need not take into 
account the values of z which satisfy simultaneously 

Z=0, Xi=0, i = l, 2, ,n. 

We also assume that there is a set of values of the variables z, x 1 , x 2 ,...., x n , 
in the vicinity of which the functions X t and Z are single-valued and analytic. 
Throughout the paper we shall confine our attention to such a domain. 

In this paper a new and complete classification is given of all the integrals 
of equation (1). In the new classification, the so-called special integrals (except 
for trivial ones) are contained in the general solution. A means is developed 
by which all these elusive special integrals can be readily determined as soon 
as the Lagrange general integral is known. 

In order to obtain the solution of (1) we take the set of equations 

dXi dz . t0 , 

-^ = -gr, t = l,2, ,n. (2) 

It is well known that (2) possesses n functionally independent integrals, one 
or more of which contain z explicitly. Let these integrals be 

^> i (z,x 1 ,x 2 , , x n )=Ci, t = l, 2, ,n, (3) 

where c x , c 2 ,...., c„ are arbitrary constants. Not all the above n equations 
are free of z. 
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Any <£,(2, x x , x 2 , . . . . , #„) =c, affords a solution of the original equation 
(1) if it involves z explicitly. Since it is an integral of the equations (2), the 
relations 



r=l CX T CZ 



i = l, 2, , n, 



(4) 



(5) 



. , <£>„, determines 



are consistent with (2). From (2) and (4) follows immediately 

,iw, X '+- S T z =°> i=1 ' 2 »• 

If the equation 

/(4>1> $2> • • • •> $>„) =0, 
in which / is an arbitrary function of the quantities',', • 
one or more values of z, then it affords solutions of (1), and is what may be 
called the Lagrange general integral. 

We have already seen that not all the solutions of every equation (1) can 
be put in this form, and since they cannot, we wish to see just how they differ 
from it. We also wish to get a solution which will be general in the true sense 
of the term. 

Let 

^{z,x ly x 2 , ,x n )=0 (6) 

be any equation whatever with the restrictions that it determines a z which 
satisfies equation (1) and that there is a set of values of the variables z, x ly 
x 2 , . . . ., x n , in the domain which we are considering, such that in their vicinity 
■^ and z are single-valued and analytic functions. 

Since ' , ', . . . ., <p n are functionally independent, not all the «-th order 
determinants of the matrix 

3' 3' 3' 3' 

dz dxi dx 2 ' ' ' ' dx n 

3' d<p 2 3<p2 d$2 

dz dx x dx 2 ' ' ' ' dx n 



3<Pn tyn ftyn <ty>n 

dz dx x dx 2 ' ' ' ' dx n 
vanish identically. 

Consequently one of the Jacobians 

j _ <? (<J>i , <?>2 , _•_•_._ ; >_$J 
<J\X l , x 2 , . . . . , x n ) 



or 



J. 



3($1, ', • • • •» 4>n) 



O ( Z, Xy , x 2 , . . . 

does not vanish identically. 
45 



^j-l > ^j + l : 



,*»)' 



5 = 1,2, 



n, 
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If J £ 0, then we can solve the system tp i =<p i (z, x 1} x 2 , . . . ., x n ), *=1, 2, 
....,«, for % , x 2 , . . . . , x n , in terms of z, ft , ft , . . . . , ft . Substituting these 
values in equation (6), we get 

^(z, x u x 2 , , x n )=g(z,$ 1 , ft, , ft) = 0. (6') 

Two cases now arise. Either z is expressible in terms of ft , ft ,...., ft , 
or equation (6') may be written in the form 

0(z)# 2 (ft, ft, , ft)=0, 

where d(z) is independent of ft, ft ,...., ft. 6(z) may be a constant. 
If z is expressible in terms of ft , ft ,...., ft , let us write 

s=0i(ft,ft,....,ft). (7) 

Differentiating z with respect to x i (»=1, 2, . . . ., n) we get 

Jl=-=it(t + t4 ^ - <") 

Multiplying the i-th equation in (8) by X { and adding the equations thus 
obtained, we get 

ix iPl = ix,i^{^ + &.,)= i i 3 M^x,+ ^X, Pl 
,=i % * % f =i i r=1 3ft\aa3 f 3e *y ^i^aftVa^ a^ ,j/t 

Hence, employing (5) and (1), we have 









or 


That is, if 


Z=0. 



> 



(9) 



2 = #l(ft, $2, > 4>») 

is a solution of equation (1), then this value of s satisfies Z—0. 

Consequently, if we wish to get solutions of form (7), we need to examine 
only those values of z which are determined by equation (9) and see whether 
they satisfy equation (1). Since no knowledge of differential equations or 
integration is involved in getting these solutions, it follows that they are trivial. 

If z is not a function of ft , ft ,...., ft , then equation (6') can be written 
in the form 

0(s)#2(ft, ft, ,ft)=0, 

where 6(z) is a function of z alone and may be a constant. 

Since the values of z determined by 6(z) =0 afford solutions of equation 
(1), we can easily determine the character of 6. The values of z determined 
by $(z) =0 are of the form0= a constant. (By exception this constant may 
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be infinity.) If z=a, a finite constant, is a solution of equation (1), a direct 
substitution in that equation shows that for s=awe have 

for all values of x 1 , x 2 , x 3 ,...., x n . 

Therefore $(z) = determines only those (non-infinite) values of z which 
satisfy Z = 0. 

If 

j __ 9 (<?>!> <?>2, ><?>») ^q 

a (z, x 1 , x 2 , . . . . , x s _ x , x s+1 , . . . . , x n ) 
where s is one of the integers 1, 2, . . . ., n, then the system <pi=fa(z, x 1 , x 2 , 
. . . ., x n ), i=l, 2, . . . ., n, can be solved for z, x x , x 2 , . . . ., x,_ x , x s+1 , . . . ., x n 
in terms of x s , fa , $ 2 , . . . . , fa . 

On substituting these values of z and x t (x s excluded) in (6), we get 

4(0, x x , x 2 , , x n )=h(x s , fa, $ 2 , , fa)=0. (6") 

Again two cases arise. Either x s is expressible in terms of fa , fa , . . . . , 

fa, or 

h(x s , fa, 4> 2 , , fa)^B s {x s )h 2 {fa, 4>2, , $«) =0, 

where 0,(« 4 ) is a function of x s alone and may be a constant. 

If x s is expressible as a function of. fa, fa, . . . . , fa, letit be written as 

a*.=M$i> $2, ,<?>«)• ( 10 ) 

Differentiating #, with respect to ^ (i=l, 2, . . . .,n), we get 

**=i!£(&+£4 '=^ - <"> 

where 8 is is equal to 1 or according as i is or is not equal to s. 

Multiplying equations (11) in order by X x , X 2 ,...., X n and adding, we get 

Hence, employing (5) and (1), we get 

dfa\ dz ~ ' dz 






or 

X s =0. (12) 

That is, if 

x s —h 1 (fa, fa, . . . .,fa) 

affords a solution of (1), then this value of x s also satisfies (12). 

Consequently we have the trivial case again. Such solutions can be deter- 
mined by processes independent of integration and involve no knowledge of 
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partial differential equations. To get these solutions we take those values of 
z which cause X s to vanish, and determine whether they satisfy equation (1). 
If x s is not a function of fa, fa, . . . ., fa, then equation (6") assumes the 
form 

0.(*.)M*i, *»,••• •>*•)=<>, 

where 6 s (x a ) is a function of x s alone and may be a constant. 

Summing up our results we have the following theorem : 

Let the general linear non-homogeneous partial differential equation be 
written in the form 

n dz 

%X i (z,x 1 ,x 2 , ,x n )p i =Z(z,x 1 ,x 2 , ,x n ), Pi=-~r- • ( A ) 

i=l OXi 

We assume that there is a set of values of the variables s, x x , x 2 ,...., x n , in 
the vicinity of which the functions Z,X X ,X 2 , . . . .,X n are single-valued and ana- 
lytic. We confine our attention to such a domain. We assume Z, X l7 X 2 , 
. . . . , X n to have no common factor. 

Solutions may be obtained by examining the values of z determined by 
Z=0, by X 1 =0, by X 2 =0, . . . ., by X n =0, and by seeing whether these satisfy 
the differential equation (A). These solutions are, however, trivial. 
All other solutions are determined as follows: 
Let n functionally independent solutions of the equations 



be 



or 



dx x dx 2 


dx n dz 


Xi X 2 


Xn ~ Z 


fa(z, x x ,x 2 , , x n )=fa, 


t=l, 2, 


One of the Jacobians 

j _ d(fa, fa, . . 
d{x 1} x z , .. 


••»*•) 

••>*») 


j . d(fa, fa, , fa) 




a (z, x x , x 2 , . . . . , x s _ x , x s+1 , . . . . 


1 «») ' 



n. 



:1, 2, ,n, 

does not vanish identically, since fa, fa, . . . ., fa are functionally independent. 

Two cases arise. 

Case I. If J £ 0, the general single-valued analytic solution, except for the 
trivial solutions afforded by 

Z=0, 
may be written in the form 

6(z)F(fa,fa,....,fa)=0, 

F being an arbitrary function of fa, fa, . . . ., fa involving z, andd(z) =0 being 
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an equation in z alone, determining only those (non-infinite) values of z which 
satisfy Z=0. 

Case II. If J s ^0 (s being an integer such that l<s<n), the general 
single-valued analytic solution, except for the trivial solutions afforded by 

x s =o, 

may be written in the form 

F being an arbitrary function 0/ $1, $ 2 » • • • -to) involving z, and Q s (x s ) being a 
function of x s alone. 

In any case the general single-valued analytic solution of (A), exclusive of 
the trivial solutions mentioned above, is 

e(z)6 1 (x 1 )6 2 (x 2 ) . . . .6 n (x n )F(to, to, ■ • ■ -, fc,)=0, 

F being an arbitrary function of to, to, ••••■> to which involves z explicitly; 
and d(z), 1 (x 1 ), B 2 (x 2 ), • • • -, B n (x n ) being functions restricted as before. 

As an illustration of the theory let us take the example 

Pi tan rCi+^2 tan #2= tan z. 
Here we get for our auxiliary solutions 

sin x 1 sin x 2 

to = — , to = — • • 

sin z T sin z 

The Lagrange general solution is then 

/(*i,*i)=0. 

where / is an arbitrary function of to and $2 involving z explicitly. 

We can easily verify that 

, . sin n+1 x 1 

V=sin z -. =0 

sin" x 2 

affords a solution of the original equation. But, since 

d(z, x l7 x 2 ) ^ ' 
4 cannot be put in the form of the Lagrange general integral. 

Since J — ~, . =£0, we know that ^ can be put in the form 

d(x x , x 2 ) ^ 

6 {e)F fa >**)•> 
thus we have 



, . / to— <K +1 \ 
^ s sm z( -^ — P — 



xl 
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We see that in this example 0(z) = sin z=0 determines only those values 
of z which cause Z, = tan z, to vanish. 
Likewise, since 

^^ and j a (ft, ft) g 

3(0, x 2 ) ^ d(z, x x ) ^ 

^ can be put in the forms 0i(fla).Fi(ft, ft) and 2 (x 2 )F 2 (<p 1 , q> 2 ) ; thus we have 

, . /I <tf\ . (1 <tf +1 \ 
^sm^---]=sm^-^]. 

An example in which 0=0 or 0^=0 may determine an infinite value of the 
variable is afforded by the equation 

x\ p 1 ~x 1 x 2 p 2 = —x\ . 

Here the auxiliary solutions are 

ft = %«2, ft = 3— 3x • 

A so-called special solution is 

^—Ssxj^ — xl=0, 

3 

= — ftft- 

x 2 

3 

Here 6 2 (x. 2 ) is- — ; equating 2 to zero, we have #2=00 . 
x 2 

An example of the trivial solutions is afforded by the equation 

(z—e x °)p 1 +p 2 =z. 

Here the auxiliary solutions are 

e Tz 
ft = — > q> 2 =x x — z + e\ 

s> 

Trivial solutions are 

^ 1= z—e x '=0, from X x =0, 
and 

^ 2 =0=O, from Z =0. 

Neither ^1 nor ^2 is a function of fy and q> 2 . 

Indiana University. 



